## Loading required package: foreach
##
## Attaching package: 'foreach'
## The following objects are masked from 'package:purrr':
##
## accumulate, when
## Loading required package: iterators
##
## Attaching package: 'data.table'
## The following objects are masked from 'package:lubridate':
##
## hour, isoweek, mday, minute, month, quarter, second, wday, week,
## yday, year
## The following objects are masked from 'package:dplyr':
##
## between, first, last
## The following object is masked from 'package:purrr':
##
## transpose
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.072208 -0.006215 -0.002112 0.005087 0.048581
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0101285 0.0002771 36.56 <2e-16 ***
## test_set$FCH4 0.2312130 0.0097955 23.60 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.009789 on 1686 degrees of freedom
## Multiple R-squared: 0.2484, Adjusted R-squared: 0.2479
## F-statistic: 557.1 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.168184 -0.004063 -0.000949 0.003350 0.090493
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0063687 0.0001099 57.96 <2e-16 ***
## result$FCH4 0.5220989 0.0040188 129.91 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.00779 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.7143, Adjusted R-squared: 0.7143
## F-statistic: 1.688e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.060437 -0.006980 -0.001889 0.004907 0.052850
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0107836 0.0002939 36.69 <2e-16 ***
## test_set$FCH4 0.2188734 0.0105320 20.78 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01049 on 1686 degrees of freedom
## Multiple R-squared: 0.2039, Adjusted R-squared: 0.2035
## F-statistic: 431.9 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.159933 -0.004229 -0.000997 0.002916 0.107940
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0066272 0.0001145 57.9 <2e-16 ***
## result$FCH4 0.5211818 0.0041861 124.5 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008114 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.6967, Adjusted R-squared: 0.6966
## F-statistic: 1.55e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.049790 -0.006012 -0.001369 0.005024 0.066616
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0097941 0.0002695 36.34 <2e-16 ***
## test_set$FCH4 0.2580954 0.0105396 24.49 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.009462 on 1686 degrees of freedom
## Multiple R-squared: 0.2624, Adjusted R-squared: 0.2619
## F-statistic: 599.7 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.124486 -0.003897 -0.000916 0.003160 0.124550
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0062906 0.0001051 59.83 <2e-16 ***
## result$FCH4 0.5428552 0.0038450 141.19 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007453 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.7471, Adjusted R-squared: 0.747
## F-statistic: 1.993e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.044136 -0.006188 -0.001447 0.004768 0.067381
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0097693 0.0002867 34.08 <2e-16 ***
## test_set$FCH4 0.2705066 0.0108225 25.00 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01004 on 1686 degrees of freedom
## Multiple R-squared: 0.2704, Adjusted R-squared: 0.2699
## F-statistic: 624.7 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.078096 -0.004011 -0.000886 0.003129 0.115947
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0062343 0.0001073 58.13 <2e-16 ***
## result$FCH4 0.5442627 0.0039225 138.75 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007603 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.7404, Adjusted R-squared: 0.7404
## F-statistic: 1.925e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.060971 -0.006431 -0.001339 0.005305 0.060552
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0097295 0.0002955 32.92 <2e-16 ***
## test_set$FCH4 0.3161960 0.0117849 26.83 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.0102 on 1686 degrees of freedom
## Multiple R-squared: 0.2992, Adjusted R-squared: 0.2988
## F-statistic: 719.9 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.080506 -0.004035 -0.000927 0.003174 0.085621
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.006229 0.000104 59.9 <2e-16 ***
## result$FCH4 0.557341 0.003803 146.5 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007372 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.7609, Adjusted R-squared: 0.7608
## F-statistic: 2.147e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.042142 -0.007056 -0.002321 0.005192 0.090158
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0107851 0.0003239 33.30 <2e-16 ***
## test_set$FCH4 0.2559779 0.0114887 22.28 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01163 on 1686 degrees of freedom
## Multiple R-squared: 0.2275, Adjusted R-squared: 0.227
## F-statistic: 496.4 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.097414 -0.004061 -0.001158 0.002898 0.132175
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0062451 0.0001204 51.87 <2e-16 ***
## result$FCH4 0.5599477 0.0044035 127.16 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008536 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.7055, Adjusted R-squared: 0.7055
## F-statistic: 1.617e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.102662 -0.006545 -0.001355 0.005129 0.046061
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0105528 0.0002852 37.0 <2e-16 ***
## test_set$FCH4 0.2198590 0.0102740 21.4 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01015 on 1686 degrees of freedom
## Multiple R-squared: 0.2136, Adjusted R-squared: 0.2131
## F-statistic: 457.9 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.225351 -0.004068 -0.000969 0.003185 0.079733
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0064394 0.0001128 57.11 <2e-16 ***
## result$FCH4 0.5270819 0.0041237 127.82 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007993 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.7077, Adjusted R-squared: 0.7076
## F-statistic: 1.634e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.071279 -0.006088 -0.001542 0.004459 0.051965
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0101719 0.0002817 36.11 <2e-16 ***
## test_set$FCH4 0.2509637 0.0105253 23.84 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.00981 on 1686 degrees of freedom
## Multiple R-squared: 0.2522, Adjusted R-squared: 0.2517
## F-statistic: 568.5 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.164506 -0.003931 -0.000896 0.003057 0.091802
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0063248 0.0001078 58.67 <2e-16 ***
## result$FCH4 0.5340847 0.0039424 135.47 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007642 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.7311, Adjusted R-squared: 0.7311
## F-statistic: 1.835e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.053052 -0.006865 -0.002164 0.004994 0.062974
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0103832 0.0002952 35.17 <2e-16 ***
## test_set$FCH4 0.2911413 0.0115847 25.13 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01029 on 1686 degrees of freedom
## Multiple R-squared: 0.2725, Adjusted R-squared: 0.2721
## F-statistic: 631.6 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.062335 -0.004168 -0.001128 0.003099 0.106795
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0063393 0.0001069 59.3 <2e-16 ***
## result$FCH4 0.5563607 0.0039096 142.3 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007578 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.75, Adjusted R-squared: 0.75
## F-statistic: 2.025e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.086313 -0.007140 -0.001803 0.005344 0.053994
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0113257 0.0002863 39.56 <2e-16 ***
## test_set$FCH4 0.1865034 0.0096318 19.36 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01038 on 1686 degrees of freedom
## Multiple R-squared: 0.1819, Adjusted R-squared: 0.1814
## F-statistic: 374.9 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.211087 -0.004384 -0.001055 0.003223 0.115594
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0068941 0.0001183 58.27 <2e-16 ***
## result$FCH4 0.4995486 0.0043271 115.45 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008387 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.6638, Adjusted R-squared: 0.6638
## F-statistic: 1.333e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.061838 -0.005811 -0.001497 0.004491 0.058570
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0099197 0.0002776 35.74 <2e-16 ***
## test_set$FCH4 0.2550037 0.0107117 23.81 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.009777 on 1686 degrees of freedom
## Multiple R-squared: 0.2516, Adjusted R-squared: 0.2511
## F-statistic: 566.7 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.136980 -0.004027 -0.000989 0.003057 0.091629
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0062858 0.0001073 58.61 <2e-16 ***
## result$FCH4 0.5448636 0.0039225 138.91 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007603 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.7409, Adjusted R-squared: 0.7408
## F-statistic: 1.929e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.048217 -0.006349 -0.001687 0.005078 0.054675
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0100743 0.0002925 34.44 <2e-16 ***
## test_set$FCH4 0.2773132 0.0115507 24.01 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01009 on 1686 degrees of freedom
## Multiple R-squared: 0.2548, Adjusted R-squared: 0.2543
## F-statistic: 576.4 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.085364 -0.003945 -0.000902 0.003066 0.093297
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0061941 0.0001067 58.07 <2e-16 ***
## result$FCH4 0.5533471 0.0039011 141.84 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007562 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.7488, Adjusted R-squared: 0.7488
## F-statistic: 2.012e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.085089 -0.006953 -0.001988 0.004998 0.071706
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0099375 0.0003258 30.5 <2e-16 ***
## test_set$FCH4 0.2958720 0.0124341 23.8 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01123 on 1686 degrees of freedom
## Multiple R-squared: 0.2514, Adjusted R-squared: 0.251
## F-statistic: 566.2 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.180234 -0.004005 -0.001030 0.003031 0.119398
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0057118 0.0001143 49.95 <2e-16 ***
## result$FCH4 0.5856925 0.0041817 140.06 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008106 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.744, Adjusted R-squared: 0.744
## F-statistic: 1.962e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.060942 -0.006939 -0.001814 0.005275 0.074804
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.010712 0.000304 35.23 <2e-16 ***
## test_set$FCH4 0.233880 0.011623 20.12 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01074 on 1686 degrees of freedom
## Multiple R-squared: 0.1937, Adjusted R-squared: 0.1932
## F-statistic: 404.9 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.133127 -0.004161 -0.000960 0.003293 0.122737
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0064488 0.0001131 56.99 <2e-16 ***
## result$FCH4 0.5418519 0.0041381 130.94 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008021 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.7176, Adjusted R-squared: 0.7175
## F-statistic: 1.715e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.067802 -0.006469 -0.001742 0.004936 0.063292
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0106747 0.0002769 38.55 <2e-16 ***
## test_set$FCH4 0.2352920 0.0101966 23.08 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.009838 on 1686 degrees of freedom
## Multiple R-squared: 0.24, Adjusted R-squared: 0.2396
## F-statistic: 532.5 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.164741 -0.004099 -0.000996 0.003139 0.123711
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0065806 0.0001095 60.08 <2e-16 ***
## result$FCH4 0.5299611 0.0040058 132.30 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007765 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.7217, Adjusted R-squared: 0.7217
## F-statistic: 1.75e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.089485 -0.006665 -0.001608 0.005269 0.056685
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.010933 0.000285 38.37 <2e-16 ***
## test_set$FCH4 0.199724 0.009266 21.55 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01038 on 1686 degrees of freedom
## Multiple R-squared: 0.216, Adjusted R-squared: 0.2156
## F-statistic: 464.6 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.205545 -0.004388 -0.001019 0.003426 0.104318
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0069286 0.0001172 59.1 <2e-16 ***
## result$FCH4 0.4906196 0.0042878 114.4 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008311 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.6599, Adjusted R-squared: 0.6598
## F-statistic: 1.309e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.061226 -0.006867 -0.001887 0.005102 0.056247
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0105388 0.0003071 34.32 <2e-16 ***
## test_set$FCH4 0.2756657 0.0118356 23.29 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01064 on 1686 degrees of freedom
## Multiple R-squared: 0.2434, Adjusted R-squared: 0.243
## F-statistic: 542.5 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.084846 -0.004159 -0.001164 0.003180 0.091598
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0063361 0.0001097 57.75 <2e-16 ***
## result$FCH4 0.5490025 0.0040127 136.82 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007778 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.735, Adjusted R-squared: 0.735
## F-statistic: 1.872e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.113537 -0.006393 -0.001591 0.005115 0.057936
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0106040 0.0002985 35.52 <2e-16 ***
## test_set$FCH4 0.2354589 0.0111378 21.14 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01047 on 1686 degrees of freedom
## Multiple R-squared: 0.2095, Adjusted R-squared: 0.2091
## F-statistic: 446.9 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.232980 -0.004157 -0.000997 0.003369 0.072216
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0064144 0.0001125 57.01 <2e-16 ***
## result$FCH4 0.5349993 0.0041151 130.01 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.007976 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.7147, Adjusted R-squared: 0.7146
## F-statistic: 1.69e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.073614 -0.006999 -0.001671 0.005542 0.066899
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.010556 0.000313 33.72 <2e-16 ***
## test_set$FCH4 0.254681 0.010930 23.30 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01111 on 1686 degrees of freedom
## Multiple R-squared: 0.2436, Adjusted R-squared: 0.2431
## F-statistic: 542.9 on 1 and 1686 DF, p-value: < 2.2e-16



## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.173372 -0.004090 -0.000957 0.003194 0.116958
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0060873 0.0001174 51.85 <2e-16 ***
## result$FCH4 0.5586618 0.0042938 130.11 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.008323 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.715, Adjusted R-squared: 0.7149
## F-statistic: 1.693e+04 on 1 and 6749 DF, p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`


## [1] "RF Summary of test data:"
##
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.103887 -0.006922 -0.002189 0.005336 0.060620
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0109562 0.0002942 37.24 <2e-16 ***
## test_set$FCH4 0.2187464 0.0100288 21.81 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.01048 on 1686 degrees of freedom
## Multiple R-squared: 0.2201, Adjusted R-squared: 0.2196
## F-statistic: 475.8 on 1 and 1686 DF, p-value: < 2.2e-16




## [1] "RF Summary of whole data:"
##
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.221616 -0.004182 -0.001106 0.003307 0.088614
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 0.0065925 0.0001145 57.56 <2e-16 ***
## result$FCH4 0.5145586 0.0041891 122.83 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.00812 on 6749 degrees of freedom
## (7936 observations deleted due to missingness)
## Multiple R-squared: 0.6909, Adjusted R-squared: 0.6909
## F-statistic: 1.509e+04 on 1 and 6749 DF, p-value: < 2.2e-16